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PROOEDE ET DISPOS1TIF D'EGAHSATION 
ET DE PECOPAGE ITERATIF 
POUR COMMUNICATIONS HAUT-PEBIT 
SUR CANAUX A ANTENNES MULTIPLES 
EN EMISSION ET EN RECEPTION 

La presente invention est relative au domaine des communications 
numeriques. 

Elle concerne plus particulierement un precede et un dispositif 
legalisation et de decodage iteratif pour communications haut-debit sur 
canaux a antennes multiples en emission et plus particulierement sur canaux 
selectifs en frequence a antennes multiples en emission et reception. 

DOMAINE TECHNIQUE GENERAL i 

Avec le developpement des technologies antennaires, les modeles de ' 
communications reposant sur le TDMA, le CDMA, I'OFDM et leurs eventuelles • 
combinaisons sont systematiquement revisites et etendus au cas de canaux ; 
MIMO. Deux grandes classes de codage espace-temps sont envisagees, qui - 
different dans leurs philosophies : 

Les transmissions en couches BLAST developpees par les Bell 
Labs ont pour but d'utiliser le multiplexage spatial pour augmenter le debit 
transmissible sur le canal. II existe deux manieres de coupler codage 
correcteur d'erreurs de canal et BLAST, conduisant a deux modeles de codage 
espace-temps en couches distinctes. Dans le premier modele (type ST BICM), 
un codage unique liant les differentes couches est applique sur les donnees 
avant entrelacement spatio-temporel et demultiplexage spatial [1, LST-I, p. 
1350]. Dans le second modele (type MAC), les donnees sont demultiplexees 
spatialement puis codees en entrelacees independamment sur chaque couche 
[1, LST-ll.p. 1350]; 
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Les schemas de diversite utiiisent la diversite spatiale pour mieux 
combattre les effets des evanouissements du canal, et ameliorer les 
performances. On distingue les STTC, qui amenent un gain de codage, et les 
STBC (et leurs generalisations OD, LDC...), qui n'amenent pas de gain de 
5 codage et qui sont souvent couples ayec les meilleures i jnpdulations codees en 
treillis connues pour le canal gaussien. 

De recentes contributions theoriques prouvent que les gains de 
multiplexage spatial et de diversite spatiale sont lies au travers d'une relation 
de compromis. Par exemple, les deux approches precitees ne permettent pas 
10 d'augmenter la capacite dans les memes proportions pour un nombre 

arbitrairement elev e d T an tennes d'e missi on M et d'antennes de reception N _On_ 

peut verifier que les schemas de diversite a remission ne sont optimaux du 
point de vue de la theorie de flnformation que dans la configuration MISO 
{M T N} = {2,1}, ce qui remet serieusement en cause leur pertinence pour les 
15 communications haut-debit sur des canaux MIMO a grand nombre d'antennes, 
Seule Tapproche ST-BICM garantit une augmentation lineaire de la 
capacite avec min {M, N}. 

C'est la raison pour laquelle il convient de privilegier ce type d'emission et 
de concentrer les efforts sur la conception d'un recepteur adapte garantissant 
20 une performance aussi proche que possible des limites fondamentales. A la 
difference des STBC, les ST-BICM n'assurent pas I'orthogonalite entre les flux 
de donnees emis par les differentes antennes. 

Dans le contexte d'une transmission sur canal MIMO selectif en 
frequence, cette strategie de codage appelle done, en reception, une detection 
25 conjointe des donnees pour supprimer a la fois interference spatiale entre 
antennes (MAI) et Interference temporelle entre symboles (ISI). 

La fonctlon de detection conjointe des donnees constitue le point delicat 
et complexe du recepteur, et ce d'autant plus que les ST-BICM necessitent, par 
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essence, un traitement iteratif entre detection et decodage de canal pour 
capturer toute la diversite spatiale. 

PRESENTATION GENERALE DE L'AR T ANTER1EUR 

Decodaae iteratif sur cananx MIMO non selectifs en fequence 
Le decodage iteratif des ST-BICM sur canaux MIMO non selectifs en 
frequence constitue un cas particulier du probleme a resoudre sous la double 
hypothese d'une elimination parfaite de I'lSI et d'une recuperation de toute 
I'energie disponible sur chaque symbole. 

Le probleme se reformule alors en ces termes : « suppression de la MAI 
sur un canal MIMO plat en presence de bruit colore spatialement ». II necessite 
le calcul des APP sur les bits des symboles des vecteurs transmis a chaque 
instant (statistique sur les donnees a transmettre au decodeur). ( 
Deux methodes ont ete retenues dans la litterature : 

Calcul exact des APP- base sur une liste exhaustive (MAP). La t 
complexite de cette methode est prohibitive, puisque exponentielle en le:. 
nombre d'antennes emettrices M et polynomiale en I'ordre Q des constellations ] 
employees. 

- Calcul approche des APP a partir d'une liste reduite de vecteurs 
candidats generee par un algorithme de decodage par sphere [2]. L'algorithme 
de decodage par sphere presente au moins deux avantages : il. preserve le 
critere MAP et possede une complexite en O (M 3 ), insensible a I'ordre Q des 
constellations employees. On precise que l'algorithme de decodage par sphere 
peut etre vu comme une variante de Talgorithme de decodage sequentiel de 
Fano. 

-• Calcul approche des APP a partir d'une liste reduite de vecteurs 
candidats generee par un algorithme sequentiel employant une pile (stack 
algorithm). 
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L'algorithme de decodage par sphere a ete intensivement etudie pour les 
canaux MIMO non selectifs en frequence. Sa generalisation au cas des canaux 
MIMO selectifs en frequence est non triviale. 

Une approche brutale serait de I'appliquer par blocs. Cette approche est 
5 peu interessante car elle augmente artificiellementjla dimension de I'espace de 
recherche (et done la complexite) et introduit un decodage redondant des 
symboles. 

L'algorithme decrit dans [2] est derive pour des canaux MIMO non 
selectifs en frequence. Sa generalisation au cas des canaux MIMO selectifs en 
io frequence est un probleme ouvert. 

La performance de ces deux algorithmes est essentiellement determinee 

par ia qualite de la liste de candidats generee, preferentiellement de taille fixe. 

II est primordial qu'elle comporte le meilleur candidat decodant un bit 1 et le 

meilleur candidat decodant un bit 0 a toute position. 
15 Decodage iteratif sur canaux MIMO selectifs en frequence 

Plusieurs types de detecteurs inseres dans des structures iteratives ont 

ete proposes. 

Detecteur a entrees et sorties ponderees realisant le critere du MAP 
base sur l'algorithme BCJR. La complexite de ce recepteur, en 0(Q UM ), est 
20 clairement prohibitive pour des modulations d'ordre eleve et des canaux MIMO 
a grand nombre d'entrees et memoire elevee. 

Detecteurs a entrees et sorties ponderees sous-optimaux bases sur 
des algorithmes de recherche efficaces dans des treillis fortement reduits en 
nombre d'etats [3], et, le cas echeant, precedes de filtres a minimum de phase. 
25 Cette approche est encore limitee par sa complexite, en O(0f- r xM ) oil L r 
represente la longueur de contrainte reduite. En particulier, I'emploi de 
constellations a grand nombre d'etats est exclu. 

Detecteurs par suppression successive ou parallele de I'interference 
[1] , utilisant des filtrages bases sur le critere MMSE (MMSE/SIC ou 
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MMSE/PIC). L'avantage principal de cette approche reside dans la faible 
complexite de mise en oeuvre. Toutefois, elle ne se revele efficace que pour 
les canaux MIMO tres peu selectifs en frequence et souffre de la mauvaise 
qualite de la premiere statistique pour les canaux MIMO fortement selectifs en 
frequence, entraTnant un demarrage tardif de ralgorithme iteratif (ecart a la 
capacite ou a la probabilite de coupure) 11 ,' Fig. 9]. 

D'autres types de detecteurs ont ete etudies isolement, i.e. en dehors 
d'un contexte de decodage iteratif. 

Detecteur iteratif par suppression d'interference . entre symboles 
applique en sortie d'un filtre adapte au canal de transmission (MF/IC) [4]. Cette 
approche, tres simple, a ete proposee dans un contexte SISO. Cependant, son 
efficacite est fortement correlee aux caracteristiques de la reponse 
impulsionnelle du canal. Pour des canaux a longues reponses impulsionnelles 
et coefficients aleatoires, le detecteur MF/IC foumit des performances proches 
du detecteur optimal MLSE. Pour des canaux a courtes reponses 
impulsionnelles et/ou taps non necessairement aleatoires, les courbes de 
performance saturent et presented des paliers a des taux d'erreurs 
relativement eleves. Ce type de phenom§ne peut egalement survenir 
accidentellement pour des reponses impulsionnelles longues, mais avec une 
probabilite faible. 

Detecteur MIMO a retour de decisions dures utilisant un filtrage 
avant calcule a partir du critere MMSE (MIMO MMSE/DFE). Cette methode 
constitue un bon compromis performance/complexite. Toutefois, I'emploi du 
DFE presente un serieux desavantage : les decisions internes ne peuvent etre 
utilisees que pour ameliorer les decisions futures compte tenu de la maniere 
sequentielle dont Palgorithme DFE traite les donnees regues. Ainsi, seute IMSI 
postcurseur peut etre soustraite au signal utile. Meme dans I'hypothese 
optimiste d'une suppression parfaite de IMSI postcurseur, la performance du 
DFE est limite par I'ISI precurseur residuel et I'augmentation du bruit. 
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PRESENTATION DE LMNVENTION 

Un but de I'invention est de proposer un recepteur avance pour les 
transmissions numeriques haut-debit pour canaux selectifs en frequence a 
5 antennes multiples en emission et reception, qui soit peu complexe - et 
notamment qui ne necessite pas des puissances de calcul importantes - tout 
en permettant de traiter Interference spatiale entre antennes et ('interference 
temporelle entre symboles. 

A cet effet, I'invention propose un dispositif d'egalisation et de decodage 

10 iteratif pour communications haut-debit sur canaux selectifs en frequence a 
antennes multiples en emission et en reception, comportant un egaliseur a 
retour de decision qui recoit les donnees des differentes antennes de 
reception, cet egaliseur comportant un filtre avant, ainsi qu'un filtre arriere 
recursif alimente par une reconstitution ponderee des donnees calculee a partir 

15 de la sortie d'un decodeur lui-meme alimente par des moyens de prise de 
decision, le dispositif comportant en outre des moyens pour soustraire la sortie 
dudit filtre arriere aux donnees en sortie du filtre avant, les donnees ainsi 
soustraites etant transmises en entree des moyens de prise de decision avec 
la sortie du decodeur, les moyens de prise de decision generant en outre une 

20 statistique transmise a destination d'un decodeur de canal a entrees et sorties 
ponderees, lesdits moyens de prise de decision prenant en compte la 
correlation spatiale du bruit en sortie des moyens de soustraction, lesdits 
moyens de prise de decision et le decodeur etant separes par des precedes 
d'entrelacement spatio-temporel au niveau binaire, 

25 caracterise en ce que le filtre avant et arriere sont adapte de facon 

iterative pour minimiser I'erreur quadratique moyenne en sortie du 
soustracteur. 
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Un tel dispositif est avantageusement complete par ies differentes 
caracteristiques prises seules ou selon toutes Jeurs combinaisons 

techniquement possibles : 

Ies moyens de decision en sortie des moyens de soustraction 
de I'egaliseur sont de type blanchiment spatial suivi d'un 
decodage par sphere ; 

Ies moyens de decision en sortie des moyens de soustraction 
de I'egaliseur sont de type SIC/PIC (suppression d'interference 
spatiale residuelle en sortie des moyens de soustraction de 
I'egaliseur, en serie et/ou en parallele) 

le blanchiment spatial est mis en ceuvre en sortie des moyens 
de soustraction de I'egaliseur 

le blanchiment spatial est mis en ceuvre par les moyens de . 
prise de decision ; 

le blanchiment spatial est mis en ceuvre par les filtres avant et , : 
arriere ; 

a partir d'une certaine iteration, le filtre avant est un filtre a. 
adapte, le filtre arriere le meme filtre adapte moins le coefficient 
central. 

L'invention conceme egalement un systeme de communication haut- 
debit sur canaux selectifs en frequence a antennes multiples en emission et en 
reception, caracterise en ce qu'il comporte en reception un dispositif 
d'egalisation et de decodage du type precite. 

Avantageusement notamment, les moyens Remission sont de type ST - 

BICM. 

En outre, l'invention propose un precede d'egalisation et de decodage 
iteratif pour communications haut-debit sur canaux selectifs en frequence a 
antennes multiples en emission et en reception, mettant en ceuvre un egaliseur 
a retour de decision qui recoit les donnees des differentes antennes de 
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reception, cet egaliseur comportant un filtre avant, ainsi qu'un filtre arriere 
recursif alimente par une reconstitution ponderee des donnees calculee a partir 
de la sortie d'un decodeur lui-meme alimente par des rnoyens de prise de 
decision, le dispositif comportant en outre des rnoyens pour soustraire la sortfe 

5 dudit filtre arriere aux donnees en sortie du filtre avant, les donnees ainsi 
soustrattes etant transmises en entree des rnoyens de prise de decision avec 
la sortie du decodeur, !es rnoyens de prise de decision generant en outre une 
statistique transmise a destination d'un decodeur de canal a entrees et sorties 
ponderees, lesdits rnoyens de prise de decision prenant en compte la 

10 correlation spatiale du bruit en sortie des rnoyens de soustraction, lesdits 
rnoyens de prise de decision et le decodeur etant separes par des precedes 
d'entrelacement spatio-temporel au niveau binaire, 

caracterise en ce que le filtre avant et arriere sont adaptes de fagon 
iterative pour minimiser Perreur quadratique moyenne en sortie du 

15 soustracteur. 



PRESENTATION DES FIGURES 

20 D'autres caracteristiques et avantages de ^invention ressortiront encore 

de la description qui suit, laquelle est purement illustrative et non limitative et 
doit etre lue en regard des dessins annexes, sur lesquels : 

- La Figure 1 illustre le concept de la technique VBLAST ; 

La Figure 2 illustre le modele general de communication connu 
25 sous le nom de ST-BICM ; 

La Figure 3 illustre Tarchitecture d'un recepteur iteratif conforme a 
un mode de realisation de invention ; 
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La Figure 4 illustre I'architecture d'un recepteur iteratif qui peut 
etre utilise a partir par exempie de ia deuxieme iteration ou des 
iterations suivantes. 

INSCRIPTION D'UM OU PHISIEURS M ODES DE REALISATION PE 
I 'INVENTION 

Rae egjs sur la techninue BLAST et le codaqe ST BICM 

On a represents sur la figure 1 le schema general d'une architecture 
BLAST. 

Les donnees codees dans un dispositif d'encodage vectoriel EV unique, 
sont transmises, apres entrelacement spatio-temporel et demuitiplexage 
spatial, .sur une pluralite d'antennes Remission TX ; a I'autre extremite des 
canaux MIMO, les signaux emis sont recus par une pluralite d'antennes de 
reception RX qui les renvoient sur des moyens de decodage DE (canaux 
MIMO) en sortie desquels on recupere les donnees. 

Dans ce qui est decrit ci-apres, on se place dans le cas de canaux MIMO 
selectifs en frequence. La presence d'interferences entre symboles augmente 
la complexite du recepteur. 

La figure 2 illustre quant a elle le modele general de communication 

connu sous le nom de ST BICM. 

A remission, les donnees sont traitees par un code correcteur de canal 1 
(codes convolutifs, turbo-codes, codes LDPC..,), puis font l'objet d'un 
entrelacement applique au niveau binaire 2, d'un demuitiplexage spatial 3, et, 
pour chaque couche, d'un precede de modulation 4-1 a 4-M. Apres un filtrage 
de mise en forme 5A a 5M, les. donnees ainsi modulees sont transmises a des 
antennes emettrices 6-1 a 6-M. 
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Description des traitements Remission et de reception 
A 1'emission, les donnees font I'objet d'un traitement ST-BICM. 
A cet effet, on met en oeuvre les etapes suivantes : 
- Reception d'un flux de donnees ^ numeriques a un debit donne.; . „ 
Application d'un codage correcteur de canal 1 generant un flux de 
donnees numeriques codees ; 

Entrelacement des donnees numeriques codees par un entrelaceur 
2 operant au niveau bit ; 

Demultiplexage spatial 3 du flux de donnees numeriques codees 
jentrelacees pour creer une pluralite .de .flux_distincts .de. donnees-numeriques 
codees appeles couches emettrices (Mflux distincts ou couches emettrices) ; 
Puis pour chacun des M flux distincts de donnees numeriques codees : 
Modulation du flux de donnees numeriques codees entrelacees 
conformement a un schema de modulation (4-1 a 4-M) pour obtenir flux de 
symboles modules ; 

Filtrage du flux module (5-1 a 5-M) ; 

Transmission du flux module sur une antenne Remission propre 
(6-1 a 6-M). 

Le recepteur du systeme de transmission est illustre sur la Figure 3. 

II comprend un egaliseur a retour de decisions defini par un filtre avant 9, 
un soustracteur 10, un algorithme prise de decisions 11 et un filtre arriere 
recursif 12. 

Le filtre arriere recursif 12 est alimente par une reconstitution ponderee 
des donnees calculee a partir de la sortie d'un decodeur 13 lui-meme alimente 
par Talgorithme de prise de decision 11. 

Le filtre avant 9 et le filtre arriere 12 sont determines de fagon iterative 
pour minimiser Terreur quadratique moyenne (MSE) de Tegaliseur MIMO, c'est- 
a-dire pour minimiser I'erreur en sortie du soustracteur 10. 
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lis utilisent a cet effet une estimation initiaie 8 du canal MIMO. 
Le traitement mis en oeuvre lest par exemple du type de celui detaille en 

annexe II. 

On notera que dans le cadre ce traitement, on travaille sur une estimation 
vectorielle de I'erreur residuel et du bruit, les filtres avant et arriere etant 
calcules par blocs pour minimiser ce vecteur. 

Ce traitement vectoriel permet un calcul simplifie. 

Les moyens de prise de decisions 11 recoivent en entree les donnees en 
sortie du soustracteur 10, ainsi que la sortie la sortie du decodeur 13. 

L'algorithme qu'ils mettent en oeuvre peut etre de different type.ll peut 
notamment etre de type SIC/PIC (suppression d'interference spatiale residuelle 
en sortie des moyens de soustraction de I'egaliseur, en serie et/ou en 
parallele) (Annexe III). 

En variante, il peut etre un algorithme de decodage par sphere. 

Un tel algorithme possede une complexite en OfM 3 ) -ouM est le nombre 
d'antennes Remission - independante de Q. 

Ceci rend possible I'utilisation de modulations a grand nombre d'etats en 

vue d'augmenter le debit. 

La sortie des moyens 11 de prise de decision est transmise a un procede 
de desentrelacement spatio-temporel 14, qui est mis en oeuvre au niveau 
binaire et qui est interpose entre lesdits moyens de prise de decision 11 et le 
decodeur 13. 

La sortie du decodeur 13 est une probability sur bits. 
Cette probability est transmise d'une part en entree de l'algorithme 11 de 
prise de decisions et d'autre part, apres entrelacement spatio-temporel 
(procede 15), a un traitement 16 de recpnstitution ponderee des donnees. 

Les donnees ponderees ainsi reconstituees sont celles qui sont 
transmises en entree du filtre aniere 12. 



( 
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La sortie du filtre arriere 12 correspond a une regeneration ponderee des 
donnees. 

Par ailleurs, on notera que Terreur correspondant a Interference 
residuelle et au bruit qui est injecte en entree du soustracteur 10, est coloree a 
5 la fois temporeliement et spatialement. 

Si la correlation temporelle a peu d'impact sur ies traitements de 
I'egaliseur, la correlation spatiale joue un role essentiel. 

C'est pourquoi, dans un mode de realisation, on realise un 
blanchissement spatial par factorisation de Cholesky. 
10 Ce blanchiment spatial est avantageusement mise en oeuvre en sortie du 
soustracteur 10, 

Dans une variante en particulier, il peut etre mise en oeuvre par 
Talgorithme 11 de decodage par sphere lui-meme. 

II peut egalement etre mis en oeuvre en etant integre aux filtres avant et 
15 arriere 9 et 12. 

Enfin, dans une variante de mise en oeuvre possible, le filtre avant 9 est 
avantageusement remplace par un filtre adapte a partir d'un certain nombre 
d'iterations donnees, par exemple a partir de la seconde iteration. 

II en est de meme du filtre arriere 12. 
20 C'est ce qu'illustre la figure 4, sur laquelle le filtre avant 9 est un filtre 

adapte MF, tandis que le filtre arriere 12 est un filtre adapte sur lequel on a 
supprime le pic correspondant au coefficient central. 
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ANNEXE I - Liste des acronymes 

TDMA : Acces multiple a repartition dans le temps (Time Division Multiple 
Access) 

5 CDMA : Acces multiple par repartition dahs le code (Code Division Multiple 
Access) 

OFDM : Muitiplexage par repartition orthogonale de la frequence (Orthogonal 
Frequency Division Multiplex) 

MIMO : Entree multiple sortie multiple (Multiple Input Multiple Output) 
10 SISO : Simple entree simple sortie (Single Input Single Output) 

BLAST : Couche spatio-temporelle des laboratoires Bell (Bell Labs Layered 
Space Time) 

ST-BICM : Modulation codee spatio-temporelle avec entrelacement binaire 
(Space-Time Bit-Interleaved Coded Modulation) 
15 STTC : Code treillis spatio-temporel (Space-Time Trellis Codes) 
STBC : Code bloc spatio-temporel (Space-Time Block Codes) 
MAI : Brouillage acces multiple (Multiple Access Interference) 
IS1 : Interference inter symbole (Inter Symbol Interference) 
APP : Probability A posteriori (A Posterior Probability) 
20 BCJR : Bahl, Cocke, Jelinek, Raviv (algorithme de) (Bahl, Cocke, Jelinek, 
Raviv (algorithm of)) 
MF : Filtre adapte (Matched Filter) 
WMF : Filtre adapte blanc (Whitening Matched Filter) 
MAP : Maximum A posteriori (Maximum A Posteriori) 
25 MMSE : Erreur quadratique moyenne minimum (Minimum Mean Square Error) 
IC : Annulation d'interference (Interference Cancellation) 
Soft-lC : Annulation d'interference souple (Soft-Interference Cancellation) 
SIC : Annulation d' interference serie (Serial Interference Cancellation) 
PIC : Annulation d' interference parallele (Parallel Interference Cancellation) 
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LDPC : Controls de parite de faible densite (Low Density Parity Check) 
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ANNEXE II 

MODEL DE COMMUNICATION 

On considere un canal P-bloc a entrees et sorties multiple a T antennes 
Remission et R antennes de reception selectif en frequence de memoire M 
(voir [2]). 

A. Modulation codee d' entrelaceme nt hinaire spatio-temporel 

Soit C un code lineaire de longueur N c et de rendementpc sur F 2 admettant un 
vecteur de donnee u e F 2 * cN et produisant un mot code c e F 2 N . Nous 
supposons que le rendement Pc inclue des bits de queue si des codes 
conventionnels sont employes. Le mot code entre dans un entrelaceur binaire 
Men concu n, dont la matrice de sortie A e F^ est segmente en P matrices 
A p e F 2 Tq * PL , p = 1 , P. Les colonnes de la matrice A p sont des vecteurs a p 

[n] e Fz Tq n = i l, que I'on nomme « vecteur d'etiquette de symbole », 

contenant T sous vecteurs aP t W 6 Ff, t = 1, .... T (un par canal d'entree), 
avec les composant binaires empiles a%>[n], ... aP<g>[n]. ou <tj> represente 
Tindex (t - 1) q+j. Dans chaque matrice A p , tous les vecteurs a p [n] sont 
modules par un modulateur sans memoire D-dimensionnel sur un ensemble de 
signal A c C D de cardinalite |A| = 2 qT a travers le codage bit a symbole u : {0, 
1} qT A . sans perte de generality nous supposons que A est un produit 
simple de T constellations complexes C de meme ordre 2* et que par 
constellation, un codage de Gray est utilise. Apres transformation du signal, ce 
processus de modulation codee peut etre aussi considere comme un schema 
de modulation codee spatio-temporelle oil tous les symboles complexes de 
mot code P*L*T sont divises en P matrices Xp e C TxL , p = 1, .... P, dont les 
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colonnes x p [n] e C T , n = 1, L sont nommees « vecteur de symbole de 
constellation ». Le « vecteur d'etiquette de symbole » a p [n] peut etre obtenu 
de x p [n] par une simple inversion de codage bit a symbole [f\ En tombant 
dans la classe generate des codes spatio-temporels, on appelle cette 
5 architecture Modulation Codee Spatio-Temporelle avec.Entrelacement.Binaire 
(STBICM) [4]. D'apres I'hypothese de filtrage ideal de Nyquist en bande limitee, 
le rendement de transmission (en bits par canal utilise) est : 

Un schema de principe relatif a ce modele de communication est montre figure 

lq 1. _ . .. 

B. Canal MIMO 

Soit H p e C R * Tx(M+1) le canal MIMO pour le bloc p, oD p = 1, P et H = {H p } 
15 I'ensemble de tous les canaux. Soient egalement X p e C TxL et Y 13 e C RxL 
respectivement le « matrice de symbole de constellation » et la « matrice de 
sortie de canal ». La sortie vectorielle de canal en bande de base equivalente a 
temps discret y 9 [n] e C R au temps n = 1 , L peut etre ecrite comme : 

y p W - {H (*)] x? [nj + [n] » 53 BP [k] \n -k] + [n] (2) 

20 Ou x? [n] e c T sont le « vecteur de symbole de constellation » transmis aux 
temps n, I'energie de chaque composante etant egale a a x 2 H p [k] e C RxT est la 
matrice du nombre de taps k de la reponse impulsionnelle du canal, w p [n] e 
C R est le vecteur de bruit complexe additif. Les vecteurs de bruit complexe 
additif w p [n] sont supposes de moyenne nulle, independant identiquement 

25 distribues, et Gaussien complexes a symetrie circulafre, et suivent done le pdf 
N (0, a 2 !). Le canal H p , constant pendant la duree du bloc correspondant, a une 
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reponse impulsionnelle finie (FIR) de longueur M+1, dont les taps de symbole 
d'espace H p [0] , . ; .,H P [M] sont les matrices aleatoires complexes R*T, de 
moyenne nutle et dont la puissance moyenne satisfait les contraintes de 
normalisation : 



:Tl' (3) 



dans le cas d'un systeme de puissance egale. L'operateur «t» correspond a 
l'operateur de transposee conjuguee. La memoire de canal egale pour tous les 
liens R*T possibles est une hypothese raisonnable etant donne que le nombre 
de composant individuel multivo.es est dicte de maniere predominante par de 
io larges structures et des objets reflechissants. 

Bloc iteratif H'annulation d'interference M1MQ 
A. Principe et ordonnancement 
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Le bloc iteratif MIMO annulateur d'interference traite chaque bloc de donnee • 

regus V\ p = 1 P . separ§ment, en utilisant Tinformation probabiliste 

aleatoire renvoyee par le decodeur de sortie . Tous les signaux et elements 
decrits sont done indexes par bloc. Nous supprimons quelque fois cette 
dependance pour simplifier les notations. Durant chaque iteration I, un filtre 
avant lineaire F 1 , applique a chaque symbole re?u de vecteur y [n], produit le 
signal de vecteur f [n]. Puis une estimation construite de maniere approprie e 
[n] du MAI et du ISI alterant x [n] est soustraite du signal de vecteur y" [n] pour 
produire 7} [n] : 

L'estimation du vecteur e' [n] pfovient de la sortie du filtre arriere B 1 excite par le 
vecteur x' [n] de tentative de decision souple sur le vecteur x[nj, etant donne la 
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15 



20 



connaissance prealable disponible (distribution de probability extrinseque) a 
Iteration I - 1 . 

Nous arrivons maintenant au coeur du turbo-egaliseur MIMO, la derivation des 
filtres a reponse impulsionnelle finie MIMO avant et arriere. 

B. Calcul des filtres avant et arriere 



Le calcul des filtres avant et arriere est inspire de [7]. Mais le contexte a 
change (contexte MIMO). Cela vaut la peine de preciser que le critere MMSE 
10 sur lequel est base cette derivation differe legerement de ceux suivis dans les 
approches concurrentes [14]. Nous commen^ons par quelques manipulations 
de base sur la sortie instantanee y^n] du filtre avant F 1 : 

Puisque : 

y M « E H W x t n ~ fe J + w M < 23 ) 

Chaque echantillon f[ri] peut etre developpe comme suit : 

^W- X) I>' [TO*] x [«-*-*] + ^ F^wfn-iJ (24) 

Une notation matricielle equivalente pour la convolution de deux filtre H c et F 1 
est : 

y* [71] - F'H^ [n] + F'w c [n] (25) 

Ou : 

F* - [PM-in] — *>] P'piftJ] €C rx ^ (26) 
est le filtre avant d'ordre L P = L Fi + U 2 + 1 , ou : 
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lie 



H[0] Hll| •-• H[M] 

H[0] H[l} -•• H [M] 

H[0] H[ll Hpfl 



€ c «I, p xTII,f+M) (27) 



H[0] H[ll— H(M) 
est la matrice de canal de bande-diagonale de Toeplitz et ou : 



x[rt] 



x Jn - ip a — M] 



(28) 



est le vecteur transmis. 
En introduisant le filtre combine : 

On obtient I'expression finale : 

; y'lrtJ^G'xcM + F'wdn] 

Avec Lgi = Lfi, Ua = L F2 + M et L 6 = L F + M. 



(30) 



(3D 



(32) 



Par anatogie, on definit le filtre arriere comme suit : 
B'= [B l [-L Bl ] 0 ...B'pfcJ] EC^ 1 " 
La sortie du correcteur de brouillage peut etre ecrite : 

2 i [n) = £ G « [k] x [n - k] | - Y, B ' W * l« " ^ + Y, * W w tn - &] . (33) 

Le vecteur de bruit et de brouillage residuel v"[n] (denomme aussi vecteur 
d'erreur de I'egaliseur MIMO) est exprime comme suit : / 
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v>] - a'[n)-G'[0]x[«] 

» ©' - *] - !>' [fcJS 1 [« - fe] + w [n-k\ 



i 



Ou 



G' - [G ? [-L Gl ] ■ - . 0 - - - G l {La,}] (35) 

5 Les filtres F 1 et B 1 sont calcules par blocs pour minimiser I'erreur quadratique 
moyenne de I'egaliseur MIMO, sous ia contrainte G'[0] = I, un probleme qui, 
dans le cas MIMO, peut etre formule de maniere compact comme suit : 

{F',BH - arg min trE<!V [n] v< fn - Is?) 
— arg min tr{Kf,| 

-[F,B}/G'[0]=I 1 J 

La minimisation MSE est executee en deux etapes successives et donne 
10 finalement : 

B ! = G' (37) 
F 1 - (E A Ht* , " 1 H e Ek)~ 1 EAl4# , " J (38) 



OCl 



etou : 

15 

avec : 



E A = [- - o 0---1 eC rxrx,c (39) 



** = [(<£ " H^Ht + a 3 !] (40) 
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<4 m 

Nous insistons sur le fait que le bruit plus le brouillage residuel a la sortie du 
soustracteur est correle spatialement et temporellement. La correlate 
temporelle n'a pas d'impact reel sur les traitements suivant, la correlation 
spatiale joue un role lei essentiel. Le bruit blanc spatial peut etre aisement 
obtenu par une simple factorisation de Cholesky de la matrice de correlat.cn 
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Proposition 2 : Puisque la matrice de correlation K' y est definie positive, la 
factorisation de Cholesky est toujours applicable. Sachant cela : 

Ou L est une matrice triangulaire inferieure. Ensuite, en appliquant F 1, = L" 1 F* et 
B '< = L; 1 B' comme filtre avant .et filtre arriere (au lieu de F» et B 1 ), la matrice de 
correlation de brouillage plus bruit egale a la matrice identite. 
Preuve : La preuve est immediate. 

K» - [(«» - 0$ L- l B*B't (L -x)t + A-¥F'» (IT 1 )*] 

• - L^^-^BW+^W'jtL-f-r'KllL- 3 )^! (43) 

Valide Talgorithme propose par la simulation de la figure 3. 

Dans la suite, nous considers toujours la sortie de I'egaliseur avec 

blanchiment spatial de sorte que : 

z l H = G z [0]x[n]-f-v l [n] 

Aveckv^letG'^L" 1 

II est interessant de remarquer que le critere MMSE correspond a une 
maximisation du SNR ^[OlX^O]}. Le SNR maximal donne pour le filtre 
adapte est atteint pour <r x 2 = • I 
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C. Decodeur APP exhaustif 

Nous definissons ici toutes les quantites probabilistes echangees au cours des 
5 iterations. La partie detection de 1'egaliseur MIMO a pour but de fournir les 
informations extrinseques sur les bits des symboies de la matrice A selon le 
nouveau modele de communication (10) ou G 1 [0] joue le role d'un canal plat 
MIMO equivalent T * T et v 1 [n] est ('interference residuelle plus bruit blanchi. 
Les rapports logarithmique APP sur tous les bits des symboies a<tj>[n] sont 
10 definis comme suit : 

M 1 J p r [o w> [n] = 0 |z* {«]] 
Par simple marginalisation, on obtient : 
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ou 



Q <S> = { a € f* «m = s } < l3 > 

Les rapports logarithmiques APP peuvent etre developpes comme 



E^w P (* l N |x=/x (a)) Pr' [a[ R ] = a] 

A 1 \n] = In * n4'^ 

M E asn ^2>(^NN=M(a))P^tafn] = a] ^ J 

20 En supposant un entrelacement spatio-temporel parfait autorise de diffuser les 
probabilites precedemment jointes dans le produit de leur composant a chiffre 
marginal : 
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p g m i* ^ w) n { t J} p-- 1 [<hm> w = «w»3 (15) 

En supposant que v' [n] a ete blanchi spatialement par un moyen quelconque 
(exemple factorisation de Cholesky), la metrique Euclidienne peut etre utilisee 
dans revaluation de la probabilite : 

p^ww^l-^w-tfW 3 *} • ° 6) 

Dans le cas general, la matrice de correlation K' v du vecteur de bruit doit etre 
prise en compte dans revaluation de la probabilite : 

p (V N |x) oc «p {- (rf W - Q l H *)* KH» (rf N - G' Pftx) } C 17) 

Comme il est etabli dans le domaine des turbo-detecteurs, les taux de 
probabilite extrinseque logarithmique sont deduit selon la formule : 

I ^>pKN 1 x ZlL (a)) n t^W> Pl ' ^ W = (18) 

Tous les echantillons d'information extrinseque logarithmique du detecteur 
MIMO sont collectes a partir de tous les blocs p = 1, .... P et rearranges, apres 
le desentrelassement spatio-temporel IT 1 , en un vecteur d'observation simple 
e R N , sur la base duquel le decodeur de sortie fourni les rapports de 
probabilite extrinseque logarithmique sur tous les bits de mots code : 



cln1 = l|£,^7{C M M>3 (19) 



Apres le desentrelacement spatSc-temporei n, le vecteur ^ c est diffuse dans P 
matrices rf* de taux de probabilite precedente logarithmique, un pour chaque 

o bloc de donnees A". Pour chaque bloc p = 1 P. les entrees de n* sont 

donnees par : < 
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4j} N = to n i r { — (20) 
Ceci permet de reecrire le vecteur x 1 [n] de prise de decision soft comme : 

*w-|E«n{i+(*« t M-i)<-'( i f*2)) m 

D. Detection M1MO a uasi-optimal par I'intermediaire du decodeur a sphere a 
liste modifiee 

La detection M1MO optimal ne peut pas etre retenue pour un scenarii de 
communication haut debit. En effet, Enumeration de point des sous- 
ensembles entiers de constellation Q <0> < fJ > et Q <1> <, i>t dont la cardinality varie 
dans 0(2 qT ), peut rapidement devenir accablante pour des modulations d'ordre 
superieur et/ou un grand nombre d'antennes de transmission. Une analyse 
soigneuse des valeurs de probabilite indique qu'un grand nombre d'entre elles 
sont negligeables. Par consequent, comme contribution significative de cet 
article, nous suggerons de remplacer Enumeration approfondie de point par 
une enumeration de point des sous-ensemble L <0> <g> et L <1> <tj> de taille 
fortement reduite, sous ensemble aussi appeles liste, qui contiennent 
uniquement les probabilites non negligeables. Les rapports de probabilite 
extrinseque logarithmique deviennent alors : 

M fn i = j ^sn<^> ( z> W |x =1* (a)) U {trJ , mtJ} IV [a ittJi) [n] = a^j 
{t ' J) n DasnW P W N F -H (a) ) 5^ IV [a (m [»] = a <ttdf> ) (18) 

En termes geometriques, contiennent des points du treillis dans une sphere de 
rayon r centree sur un point bien choisi (par exemple le point ML sans 
coritrainte ou le point ML lui-meme). Des versions modifiees de decodeur a 
sphere ont ete employees avec un certain succes pour exhiber ces listes [15]. 
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Le choix du rayon de sphere r pilote la qualite et la complexite du detecteur 
MIMO correspondant. 

F. Algorithms de detection PIC 

Le decodeur a sphere peut egalement etre considere comme trop complexe. 
Une maniere possible d'alleger encore la complexite du detecteur est de 
generer I'APP log independamment par dimension. La variable de decision 
pour le composant x est, a Iteration I : 
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Ou g k est la colonne k de la matrice G1[0]. Les rapports APP log £<y>[nl, j- 
1 ..... q sont alors calcules en supposant que r\[n] est une variable gaussienne 
de moyenne x tl et de variance estimee : 

Le premier detecteur ^iteration peut etre ameliore en utilisant une technique 
SIC de la facon suivante : 
• Initialisation 

Classifier les composants par ordre decroissant en fonction de leur 
rapport signal sur parasite (SIR), par exemple : 

, SIR{xi)>SIR(x2)>'->SXR(x T ) 
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gfgi 



SMI*) - f ' 



Recursivite 
Pourt=1 a T 
Calcul de 



Calcul du rapport APP log £'<tj>[n] , j= 1 q en supposant que r^n] est une 
variable gaussienne de moyenne x t , et de variance estim.ee : 

Calcul Pestimation MMSE de x* connaissant r^n] : 

£t[n] ~ E {x t {n]\rl[n]} 

Fin pour 
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REVINDICATIONS 

1 . Dispositif d'egalisation et de decodage iteratif pour communications haut- 
debit sur canaux selectifs en frequence a antennes multiples en emission 
et en reception, comportant un egaliseur a retour de decision qui regoit les 
donnees des differentes antennes de reception, cet egaliseur comportant 
un filtre avant (9), ainsi qu'un filtre arriere (12) recursif alimente par une 
reconstitution ponderee des donnees calcuiee a partir de la sortie d'un 
decodeur (13) lui-meme alimente par des moyens de prise de decision 
(11), le dispositif comportant en outre des moyens pour soustraire (10) la 
sortie dudit filtre arriere (12) aux donnees en sortie du filtre avant (9), les 
donnees ainsi soustraites etant transmises en entree des moyens de prise 
de decision (11) avec la sortie du decodeur (13), les moyens de prise de 
decision (1 1 ) generant en outre une statistique transmise a destination d'un 
decodeur de canal a entrees et sorties ponderees, lesdits moyens de prise 
de decision (11) prenant en compte la correlation spatiale du bruit en sortie 
des moyens de soustraction (10), lesdits moyens de prise de decision (11) 
et le decodeur (13) etant separes par des precedes d'entrelacement 
spatio-temporel au niveau binaire, 

caracterise en ce que les filtres avant (9) et arriere (12) sont adaptes de 
fa?on iterative pour minimiser I'erreur quadratique moyenne en sortie du 
soustracteur (10). 

2. Dispositif selon la revendication 1, caracterise en ce que les moyens de 
decision (11) en sortie des moyens de soustraction (10) de I'egaliseur sont 
de type blanchiment spatial suivi d'un decodage par sphere ; 

3. Dispositif selon la revendication 1, caracterise en ce que les moyens de 
decision (11) en sortie des moyens de soustraction (10) de I'egaliseur sont 
de type a suppression ^interference spatiale residuelle en sortie des 
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moyens de soustraction (10) de I'egaliseur, en serie et/ou en 
parallele(SIC/PlC) 

4. Dispositif selon Tune des revendications precedentes, caracterise en ce que 
le blanchiment spatial est mis en ceuvre en sortie des moyens de 
soustraction (10) de I'egaliseur 

5. Dispositif selon la revendication 4, caracterise en ce que le blanchiment 
spatial est mis en ceuvre par les moyens de prise de decision (11). 

6. Dispositif selon la revendication 2, caracterise en ce que le blanchiment 
spatial est mis en ceuvre par les filtres avant (9) et arriere (12); 

7. Dispositif selon l'une des revendications precedentes, caracterise en ce que 
a partir d'une certaine iteration, le filtre avant (9) est un filtre adapte, le filtre 
arriere (1 2) le meme filtre adapte moins le coefficient central. 

8. Systeme de communication haut-debit sur canaux selectifs en frequence a 
antennes multiples en 6mission et en reception, caracterise en ce qu'il 
comporte en reception un dispositif d'egalisation et de decodage selon 
l'une des revendications precedentes. 

9. Systeme selon la revendication 7, caracterise en ce que les moyens 
d'emission sont de type ST -BICM. 

10. Precede d'egalisation et de decodage iteratif pour communications haut- 
debit sur canaux selectifs en frequence a antennes multiples en emission 
et en reception, mettant en ceuvre un egaliseur a retour de decision qui 
recoit les donnees des differentes antennes de reception, cet egaliseur 
comportant un filtre avant (9), ainsi qu'un filtre arriere (12) recursif alimente 
par une reconstitution ponderee des donnees calculee a partir de la sortie 
d'un decodeur (13) lui-meme alimente par . des moyens de prise de 
decision (11), le dispositif comportant en outre des moyens (10) pour 
soustraire la sortie dudit filtre arriere (12) aux donnees en sortie du filtre 
avant (9), les donnees ainsi soustraites etant transmises en entree des 
moyens de prise de decision (11) avec la sortie du decodeur (13), les 
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moyens de prise de decision (11) generant en outre une statistique 
transmise a destination d'un decodeur de canal a entrees et sorties 
ponderees, lesdits jnoyens de prise.de decision (11) prenant en compte la 
correlation spatiale du bruit en sortie des moyens de soustraction (10), 
5 lesdits moyens de prise de decision (11) et le decodeur (13) etant separes 

par des precedes d'entrelacement spatio-temporel au niveau binaire, 
caracterise en ce que le filtre avant (9) et arriere (12) sont adaptes de 
fagon iterative pour minimiser I'erreur quadratique moyenne en sortie du 
soustracteur (10). 

10 11. Precede selon la revendication 10, caracterise en ce que le filtrage recursif 
de legalisation a retour de decisions est egalement adapte de fagon 
iterative pour minimiser I'erreur quadratique moyenne en entree de 
I'egaliseur. 

IZProcede selon Tune des revendications 9 ou 10, caracterise en ce que Ton 
15 met en ceuvre un blanchiment spatial sur les donnees soustraites 

13. Precede selon Tune des revendications 9 a 11, caracterise en ce que 
Palgorithme de prise de decision (11) est de type decodeur par sphere. 

14. Procede selon les revendications 1 1 et 12 prises en combinaison, 
caracterise en ce que le blanchiment spatial est mis en oeuvre par 

20 l'algorithme de prise de decision (1 1 ). 

15. Procede selon Tune des revendications precedentes, caracterise en ce qu'a 
partir d'une certaine iteration, le filtre avant (9) est un filtre adapte. 
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